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Abstract 



We establish the expressibihty in fixed-point logic with counting (FPC) 
of a number of natural polynomial-time problems. In particular, we show 

I I that the size of a maximum matching in a graph is definable in FPC. 

^^ This settles an open problem first posed by Blass, Gurevich and She- 

^ lah [BGS99| . who asked whether the existence of perfect matchings in 

general graphs could be determined in the more powerful formalism of 
choiceless polynomial time with counting. Our result is established by 
showing that the ellipsoid method for solving linear programs can be im- 
plemented in FPC. This allows us to prove that linear programs can 

^~~^ be optimised in FPC if the corresponding separation oracle problem can 

'^ be defined in FPC. On the way to defining a suitable separation oracle 

■^^ for the maximum matching problem, we provide FPC formulas defining 

Ji_ maximum flows and canonical minimum cuts in capacitated graphs. 



*Research supported by EPSRC grant EP/H026835. An extended abstract of this paper 
\] will appear in the proceedings of LICS 2013. 
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1 Introduction 

The question of whether there is a logical characterisation of the class P of prob- 
lems solvable in polynomial time, first posed by Chandra and Harel |CH82| . has 
been a central research question in descriptive complexity for three decades. At 
one time it was conjectured that FPC, the extension of inflationary fixed-point 
logic by counting terms, would suffice to express all polynomial-time properties, 
but this was refuted by Cai, Fiirer and Immerman |CFI92j . Since then, a num- 
ber of logics have been proposed whose expressive power is strictly greater than 
that of FPC but still contained within P. Among these are FPR, fixed-point 
logic with rank operators JDGHL09] . and CPT(Card), choiceless polynomial 
time with counting |BGS99[ IBGS02| . For both of these it remains open whether 
their expressive power is strictly weaker than P. 

Although it is known that FPC does not express all polynomial-time com- 
putable properties, the descriptive power of FPC still forms a natural class 
within P. For instance, it has been shown that FPC can express all polynomial- 
time properties on many natural graph classes, such as any class of proper 
minor-closed graphs |GrolO| . Delimiting the expressive power of FPC there- 
fore remains an interesting challenge. In particular, it is of interest to establish 
what non-trivial polynomial-time algorithmic techniques can be expressed in 
this logic. The conjecture that FPC captures P was based on the intuition that 
the logic can define all "obvious" polynomial-time algorithms. The result of Cai 
et al. and the subsequent work of Atserias et al. [ABDOQ^ showed that one impor- 
tant technique — that of Gaussian elimination for matrices over finite fields — is 
not captured by FPC. The question remains what other natural problems for 
which membership in P is established by non-trivial algorithmic methods might 
be expressible in FPC. 

For instance, it was shown by Blass et al. |BGS02| that there is a sentence 
of FPC that is true in a bipartite graph G if, and only if, G contains a perfect 
matching. They posed as an open question whether the existence of a perfect 
matching on general graphs can be defined in CPT(Card) (see also |BG05[ 
IRoslOj for more on this open question). Indeed, this question first appears 
in [BGS99J where it is stated that it seems "unlikely" that this problem can be 
decided in CPT(Card). One of our main contributions in this paper is to settle 
this question by showing that the size of a maximum matching in a general 
graph can be defined in FPC (and therefore also in CPT(Card)). 

On the way to establishing this result, we show that a number of other 
interesting algorithmic problems can also be defined in FPC. To begin with, 
we study the logical definability of linear programming problems. Here we show 
that there is a formula of FPC which defines on a polytope a point inside the 
polytope which maximises a given linear objective function, if such a point 
exists. Here, by a polytope we mean a convex set in Euclidean space given by 
finite intersections of linear inequalities (or constraints) over a set of variables, 
suitably represented as a relational structure without an ordering on the sets of 
variables or constraints. 

More specifically, we consider representations where, as in many applications 
of linear programming, the set of constraints is not given explicitly (indeed, it 
may be exponentially large) but is determined instead by a separation oracle. 
This is a procedure which, given a candidate point x, determines whether x is 
feasible and, if it is not, returns a constraint that is violated by x. It is well 



known that Khachiyan's polynomial-time algorithm for linear programming — 
the ellipsoid method — can be extended to prove that the linear programming 
and separation problems are polynomial-time equivalent (c.f., |Kha8 . GLS81] 
IGLS88| ). We show an analogous result for FPC: if a separation oracle for a 
polytope is expressible in FPC, then linear programming on that polytope is 
definable in FPC. Informally speaking, the idea is the following. Although the 
set of variables is not inherently ordered, the separation oracle induces a natural 
equivalence relation on these variables whereby two variables are equivalent if 
they cannot be distinguished in any invocation of the oracle. Given a linear 
ordering on these equivalence classes, we can define in FPC a reduction of 
the optimisation problem to an instance with an ordered set of variables by 
taking the quotient of the polytope under the induced equivalence relation. We 
show that solving the optimisation problem on this quotiented polytope — now 
using the classical polynomial-time reduction made possible by the ordering — 
allows us to recover a solution to the original problem. In practice, neither the 
equivalence classes nor the order are given beforehand; rather they are iteratively 
refined via the invocations of the separation oracle made while optimising over 
the quotiented polytope. The details of this result are presented in Section |4] 

Thus to express a problem modelled by a linear program in FPC it suffices 
to express a separation oracle in FPC. A key difficulty to expressing separation 
oracles in FPC is that a particular violated constraint must be chosen. We show 
in Section [3] that when the constraints are given explicitly, a canonical violated 
constraint can be defined by taking the sum of all the violated constraints. This 
implies that the class of feasible linear programs, when explicitly given, can be 
expressed in FPC. When the constraints are not given explicitly, it may still be 
possible to express canonical violated constraints (and hence separation oracles) 
by using domain knowledge. 

As a first application of the FPC-definability of explicitly-given linear pro- 
grams, we show that a maximum flow in a capacitated graph is definable in 
FPC. Indeed, this follows rather directly from the first result, since the flow 
polytope is of size polynomial in G and explicitly given, and hence a separation 
oracle can be easily defined from G in FPC. These results are presented in 
Section [Sj 

Next, we use the definability of maximum flows to show that minimum cuts 
are also definable in FPC. That is, in the vocabulary of capacitated graphs, 
there is a formula which defines a set of vertices G corresponding to a minimum 
value cut separating s from t. The cut G defined in this way is canonical in 
a strong sense, in that we show that it is the smallest (under set-inclusion) 
minimum cut separating s from t. The definition of minimum cut and useful 
variants are presented in Section [61 

Finally, we turn to the maximum matching problem. For a graph G = (V, E), 
the matching polytope is given by a set of constraints of size exponential in the 
size of G. We show that there is a separation oracle for this set definable from 
G in FPC, using the definability of minimum cuts. To be precise, we use the 
fact that a separation oracle for the matching polytope can be obtained from 
a computation of minimum odd-size cuts in a graph |PR82j . In Section m we 
prove that there is always a pair of vertices s, t such that a canonical minimum 
(s,i)-cut is a minimum odd-size cut. This, combined with the definability of 
canonical minimum cuts, gives us the separation oracle for matching that we 
seek. Note that it is not possible in general to actually define a canonical 



maximum matching. To see this, consider Kn, the complete graph on n vertices. 
This graph contains an exponential number of maximum matchings and for any 
two of these matchings, there is an automorphism of the graph taking one to 
the other. Thus, it is not possible for any formula of FPC (which is necessarily 
invariant under isomorphisms) to pick out a particular matching. What we can 
do, however, is to define a formula that gives the size of the maximum matching 
in a graph. This, in turn, enables us to write a sentence of FPC that is true in a 
graph G if, and only if, it contains a perfect matching. Our results on matchings 
are presented in Section [8| 

2 Background 

We write [n] to denote the set of positive integers {0, . . . , n — 1}. Given sets / 
and A, a column vector u over A indexed by / is a function u : I -^ A, and 
we write A for the set of all such vectors. Similarly, an /, J-matrix over A is 
a function M : I x J ^ A and we write Mij for M{i,j) and M^ for the row 
(vector) of M indexed by i. For an integer z, \z\ denotes its absolute value. For 
a vector v S Q^, Halloo '■— maxjg/ \vi\ denotes its infinity norm. 

2.1 Logics and Structures 

A relational vocabulary r is a finite sequence of relation and constant symbols 
(i?i, . . . , Rk, ci, . . . , q), where every relation symbol Ri has a fixed arity a^ £ N. 
A structure A = (dom(A), Rf-, . . . , R^, cf-,. . . ,cf-) over the vocabulary r (or 
a T-structure) consists of a non-empty set dom(A), called the universe of A, 
together with relations Rf^ C dom(A)°' and constants c^ S dom(A) for each 
\ < i < k and 1 < j < H. Members of the set dom(A) are called the elements 
of A and we define the size of A to be the cardinality of its universe. In what 
follows, we often consider multi-sorted structures. That is, dom(A) is given 
as the disjoint union of a number of different sorts. In this paper we consider 
only finite structures, that is structures over a finite universe. For a particular 
vocabulary r we use fin[r] to denote the set of all finite r-structures. 

Fixed-point logic with counting. Fixed-point logic with counting (FPC) 
is an extension of inflationary fixed-point logic with the ability to express the 
cardinality of definable sets. The logic has two types of first-order variable: 
element variables, which range over elements of the structure on which a formula 
is interpreted in the usual way, and number variables, which range over some 
initial segment of the natural numbers. We traditionally write element variables 
with lower-case Latin letters x,y, . . . and use lower-case Greek letters /i, 77, . . . 
to denote number variables. 

The atomic formulas of FPC[r] are all formulas of the form: /i = 77 or 
fi < rj, where /z, rj are number variables; s = t where s, t are element variables or 
constant symbols from t; and R{ti, . . . ,tm), where each ti is either an element 
variable or a constant symbol and _R is a relation symbol of arity m. The set 
FPC[t] of FPC formulas over r is built up from the atomic formulas by applying 
an inflationary fixed-point operator [ifp^^ -(/)] (t) ; forming counting terms #x'/>j 
where </> is a formula and x an element variable; forming formulas of the kind 
s = t and s < t where s, t are number variables or counting terms; as well as the 



standard first-order operations of negation, conjunction, disjunction, universal 
and existential quantification. Collectively, we refer to element variables and 
constant symbols as element terms, and to number variables and counting terms 
as number terms. 

For the semantics, number terms take values in [n + 1] and element terms 
take values in dom(A) where n := |dom(A)|. The semantics of atomic formulas, 
fixed-points and first-order operations are defined as usual (c.f., e.g., |EF99j for 
details), with comparison of number terms /i < ?7 interpreted by comparing 
the corresponding integers in [n -|- 1] . Finally, consider a counting term of the 
form 4j=x4>^ where </) is a formula and x an element variable. Here the intended 
semantics is that 4J^x4> denotes the number (i.e., the element of ['t.+I]) of elements 
that satisfy the formula (p. 

In general, a formula (f){x, fl) of FPC defines a relation over dom( A) l+l [rt + 1] 
that is invariant under automorphisms of A. For a more detailed definition of 
FPC, we refer the reader to [EF^lLibOi] . 

It is known, by the results of Immerman and Vardi |Imm861 IVar82j , that 
every polynomial-time decidable property of ordered structures is definable in 
fixed-point logic, and therefore also in FPC. Here, an ordered structure is 
one which includes a binary relation which is a linear order of its universe. 
Throughout this paper, we refer to this result as the Immerman- Vardi theorem. 

Logical interpretations. We frequently consider ways of defining one struc- 
ture within another in some logic L, such as first-order logic or fixed-point logic 
with counting. Consider two vocabularies a and r and a logic L. An m-ary 
Ij-interpretation of t in a is a sequence of formulae of L in vocabulary a con- 
sisting of: (i) a formula S{x); (ii) a formula e{x, y); (iii) for each relation symbol 
R € T oi arity fc, a formula (j)ji{xi, . . . , Xfc); and (iv) for each constant symbol 
c G T, a formula 7c(af), where each x, y or Xi is an ?7i-tuple of free variables. 
We call TO the width of the interpretation. We say that an interpretation Q 
associates a r-structure B to a cr-structure A if there is a surjective map h from 
the TO-tuples {a G (dom(A) l+l [n -f 1])™ | A |= 6[a\} to B such that: 

• h{ai) = ^1(02) if, and only if, A |= £[01,(12]; 

• i?^(/i(ai), . . . , h{ak)) if, and only if, A ^ 4>R[ai,. . . , Ofe]; 

• h{a) = c^ if, and only if, A |= 7c [a]- 

Note that an interpretation Q associates a r-structure with A only if e defines an 
equivalence relation on (dom(A) l±l [n -|- 1])™ which is a congruence with respect 
to the relations defined by the formulae (j)^ and 7c. In such cases, however, B 
is uniquely defined up to isomorphism and we write 0(A) := B. 

It is not difficult to show that formulas of FPC compose with reductions in 
the sense that, given an interpretation of cr in t and a u-formula (/>, we can 
define a t- formula 0' such that A\= (j)' if, and only if, 0(A) |= </> (see |Imm991 
Sec. 3.2]) In particular, if 6(A) is an ordered structure, for all A, then by the 
Immerman- Vardi theorem above, for any polynomial-time decidable class C, 
there is an FPC formula 4> such that A |= if, and only if, 0(A) S C. 



2.2 Numbers, Vectors and Matrices 

Let z be an integer, h > [log2(|2;|)] , B = [b] and write hit{x,k) to denote 
the k-th least-significant bit in the binary expansion of x e N. We view the 
integer 2 = s • x as a product of a sign s G {—1, 1} and a natural number 
X. We can represent z as a single-sorted structure B on a domain of bits 
B over the vocabulary tj, := {X, S,<b}- Here <b is interpreted as a linear 
ordering of B, the unary relation S indicates that the sign s of the integer 
is 1 if S^ = and —1 otherwise, and the unary relation X is interpreted as 
X^ ^ {k€ B \ bit(a;, k) = 1}. That is k e X^ when the "the fc-th bit in the 
binary expansion of x is 1." Similarly we consider a rational number q ~ s ■ ^ 
as a structure on the domain of bits B over tq := {X, D, S,<b}, where X and 
S are as before and D is interpreted as the binary encoding of the denominator 
d when D^ ^ 0. 

We now generalise these notions and consider unordered tensors over the 
rationals (the case of integers is completely analogous) . Let Ji , . . . , J^ be a 
family of finite non-empty sets. An unordered tensor T over Q is a function 
T : Ji X • • • X Jr — > Q. We write iji...j,, — gji...j,- d^ ^' to denote the element of 
T indexed by (j'l, . . . , j^) G J\ x • • • x J^. Writing to € N for the the maximum 
absolute value of integers appearing as either numerators or denominators of 
elements in the range of T, let b > [log2(|m|)] and B = [b]. The tensor T 
is then an (r -|- l)-sorted structure T with r index sorts Ji, . . . , J^ and a bit 
sort B over the vocabulary atcn.r '■= {X, D, S,<b}- Here <b is interpreted as 
before, the (r -I- l)-ary relation S is interpreted as indicating the value of the 
sign Sj^,,,j^ £ {^Ij 1} a-s before, the (r + l)-ary relation X is interpreted as 

{{ji,...,jr,k) e JiX ■■■ X Jr X B\ hit{Xj^. ..j^,k) = 1}, 

and the (r+l)-ary relation D is similarly interpreted as the binary representation 
of the denominators of T. We are only interested in the case of rational vectors 
and matrices and so define the vocabularies Tvoc := crten,! and Tmat := crten,2- 

In |HollO| it is shown that a variety of basic linear-algebraic operations 
on rational vectors and matrices described in this way can be expressed in 
fixed-point logic with counting. These include computing equality, norms, dot 
product, matrix product, determinant and inverse. 

2.3 Linear Programming 

We recall some basic definitions from combinatorics and linear optimisation. For 
further background, see, for example, the textbook by Grotschel et al. |GLS88| . 

Polytopes. Consider the rational Euclidean space Q^ indexed by a set V. 
The solutions to a system of linear equalities and inequalities over Q^ is the 
intersection of some number of half-spaces of the kind {x G Q^ | a^ x < b} spec- 
ified by the constraint a^ x < b, where a G Q^ and b £ Q. A (rational) polytope 
is a convex set P C Q^ which is the intersection of a finite number of half- 
spaces. That is to say, there are a set of constraints C, a constraint matrix 
A G Q^""^ and vector b S Q'^, such that P = PA,b ■= {x £ Q^ \ Ax < b}. 

Polytopes have an alternative characterisation as a combination of convex 
hulls and cones. Let S* be a finite set of points in Q^ and define the convex hull 



of S* 
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If P is a polytope in Q^, then there exist finite sets iS'i,S'2 C Q^ such that 
P = Psi,S2 := conv(S'i) +cone(5'2) = {xi + X2 | xi S conv(5'i),X2 G cone(S'2)}. 
The size, or bit complexity, of a vector c € Q^ (denoted (c) ) is the number 
of bits required to encode the components of c in some standard encoding of 
rational numbers. Note that (c) is at least \V\. The size of a constraint a^x < b 
is then ((^))- If Ax < & is a system of linear inequalities then its size is the 
maximum over the sizes of its individual constraints. Note that this measure 
is explicitly independent of the number of constraints in the system. The facet 
complexity {P)f of a polytope P is the minimum over the sizes of the systems 
Ax < b such that P — PA,b- The vertex complexity {P)v of a polytope P is the 
minimum over the maximum size vector in the union of sets Si, S2 Q Q^ such 
that P — conv(S'i) + conc(S'2). The facet and vertex complexity of a polytope 
are closely related: (P)v < 4|F|2(P)f and (P)f < 3\V\^{P)^ |GLS88l Lemma 
6.2.4]. 

Problems on polytopes. We are interested in two main combinatorial prob- 
lems on polytopes: linear optimisation and separation. 

Problem 1 (Linear Optimisation). Let V be a set, P C Q^ be a polytope and 
c E Q^. The linear optimisation problem on P is the problem, of determining 
either (i) an element y E P such that c^y = maxjc^a; | x £ P}, (ii) that P = 
or (Hi) that P is unbounded in the direction of c. 

An instance of the linear optimisation problem is called a linear program 
and the linear function x 1—^ c^ x is called the objective function. Over the 
years, a number of algorithms for solving linear programs have been studied. 
Early work by Dantzig |Dan63| gave a combinatorial algorithm — the simplex 
method — which traverses the vertices (extremal points) of the polytope favour- 
ing vertices that improve the objective value. Although the simplex method is 
useful in practice, it tends not to be theoretically useful because strong worst- 
case performance guarantees are not knowrFlAseries of works studying linear 
programming from a geometric perspective |Sho721 IYN76[ ISho77| culminated in 
the breakthrough of Khachiyan jKha79irKha80) which established a polynomial- 
time algorithm — the ellipsoid method — for solving linear programs. One of the 
strengths of the ellipsoid method is that it can be applied to linear programs 
where the constraints are not given explicitly. In such implicitly-defined lin- 
ear programs, we are instead given a polynomial-time algorithm, known as a 
separation oracle, for solving the following "separation problem" . 



^ Stronger guarantees are known for the average-case and smoothed complexity of the sim- 
plex method |ST04| . 



Problem 2 (Separation). Let V be a set, P C Q^ be a polytope and y G Q^. 
The separation problem on P is the problem of determining either (i) that y (z P 
or (a) a vector c € Q^ with c^y > maxjc^o; | x € P} and \\c\\oo — 1- 

Over families of rational polytopes, the optimisation and separation prob- 
lems are polynomial-time equivalent (c.f., e.g., |GLS88[ Theorem 6.4.9]). Here 
the time bound is measured in the size of the polytope and all other parameters 
of the problem. 

2.4 Representation 

When we deal with polytopes as objects in a computation, we need to choose 
a representation which gives a finite description of a polytope. In particular, in 
dealing with logical definability of problems on polytopes, we need to choose a 
representation of polytopes by relational structures. 

Definition 3. A representation of a class V of polytopes is a relational vo- 
cabulary T along with an onto function v : fin[T] — >■ V which is isomorphism 
invariant, that is, A = B implies i'(A) = i^(B). 

For concreteness, consider the vocabulary r :— Tmat W ''vec obtained by tak- 
ing the disjoint union of the vocabularies for rational matrices and vectors. A 
T-structure over a universe consisting of a set V of variables and a set C of 
constraints describes a constraint matrix A £ Q'-'^'^ and bound vector b £ Q'^' . 
Thus, the function taking such a structure to the polytope PA,b is a representa- 
tion of the class of rational polytopes. We call this the explicit representation. 

Note that the explicit representation of polytopes has the property that both 
the size of the polytope (i.e., the maximum size of any constraint) and the num- 
ber of constraints of ^{A.) are polynomially bounded in the size of A. We will 
also be interested in representations v where the number of constraints in ^{A.) 
is exponential in |A|, but we always confine ourselves to representations where 
the size of the constraints is bounded by a polynomial in A. We formalise this 
by saying that a representation v is well described if there is a polynomial p such 
that {i/{A)) — p{\A\), for all r-structures A. In particular, in all representations 
we consider the dimension of the polytope j^(A) is bounded by a polynomial in 
|A|. 

We are now ready to define what it means to express the linear optimisation 
and separation problems in FPC. 

Definition 4. We say that the linear optimisation problem for a class of poly- 
topes V is expressible in FPC with respect to a representation v : fin[T] -^ V if 
there is an FPC interpretation ofTQ^T^ec 'in T^Ty^c which takes a r-structure A 
and a vector c to a rational f and vector y such that either (i) / = 1, and ^{A) 
is unbounded in the direction of c, or (ii) f — 0, and ^{A) 7^ iffy'E v{A) and 
c^y = maxjc^a; | x e ^{A)}. 

Definition 5. The separation problem for a class of polytopes V is expressible 
in FPC with respect to a representation v : fin[T] -^ V if there is an FPC 
interpretation of Ty^c ii^ T^Ty^c which takes a structure coding a T-structure A 
and a vector y to a vector c such that either (i) y G '^(A) and c = 0, or (ii) 
c £ Q^ with c^y > maxjc^a; | x £ i^(A)} and \\c\\^ — 1. 



Figure 1: A separation oracle for explicitly- represented rational polytopes. 
A(A,b,x) 

Input: A e q^''^, b€Q^ andxe Q^ . 

Output: c G Q^ solving the separation problem for the polytope P^.t and x. 

1: if Ax < b then return 0^. 

2: Select k G C such that Aj^x > h^. 



3 Expressing the Separation Problem in FPC 

Let A e Q*^^^ be a constraint matrix and b e Q^ a constraint vector of the 
polytope Pam- Figure [I] presents a straightforward algorithm (A) for solving 
the separation problem for the explicitly-represented polytope PA,b- It is not 
hard to see that the algorithm A can be implemented in time polynomial in the 
size of the explicit natural representation of inputs A, b and x. 

If we try to express the algorithm A in fixed-point logic with counting, we 
first note that we can define in FPC all the relevant manipulations on rational 
values, vectors and matrices, such as norms, addition and multiplication |HollO| . 
even when they are indexed by unordered sets. This shows that both linesfTland 
[3] of the algorithm can be simulated in FPC. However, line [2] poses a problem as 
the logic is in general not able to choose a particular element from an unordered 
set. Our key observation here is that linearity implies that the sum of all such 
violated constraints is itself a violated constraint for non-empty polytopes and 
hence the choice made by A is superfluous. This can be formally stated as 
follows. 

Proposition 6. Let A € Q<^^^, b € Q^', a; e Q^ and C 3 S" 7^ 0. Suppose 
PA.b is non-empty and {Ax)s ^ bs for all s £ S. Define as ■— J2ses^^- Then 
OgX > maxjajy | y G PA,b} and as y^ 0^ ■ 

Proof. Define bs '■— "^ses^s- That ajx > bs is immediate from linearity. 
Since the polytope is non-empty pick any point y £ Pa, 6- By definition. Ay < b. 
Linearity implies that ajy < 65. Thus ajx > 65 > maxjajy | y £ PA.b}- This 
also implies that as y^ 0^ . D 

This observation leads to a definition in FPC of the separation problem for 
PA.b with respect to x. Specifically, let 5' C C be the set of constraints which 
violate the inequality Ax < b. This set can be defined by a FPC formula using 
rational arithmetic. If S is empty, expressing c = 0^ correctly indicates that 
X G PA,b- Otherwise S is non-empty; let as be the sum of the constraints which x 
violates. Since the set 5* is definable in FPC so is the sum of constraints indexed 
by S. If as 7^ 0^, Proposition [6] implies that expressing c as the division of as 
by its (non-zero) infinity norm correctly indicates a separating hyperplane for 
PA.b through x; moreover, both operations are in FPC. Otherwise, as = 0^ and 
Propositionl6]indicates that PA.b is empty. This means that any non-zero vector 
defines a separating hyperplane for Pa, 6- Thus it suffices for the interpretation 
to express the vector c — 1^ . Overall, the above discussion gives us a proof of 
the following theorem. 
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Theorem 7. There is an FPC interpretation of Ty^c in Tmat W T^ec W Ty^c ex- 
pressing the separation problem for the class of polytopes explicitly given by 
constraints and represented naturally as Tmat W T^ec-structures. 

4 Reducing Optimisation to Separation in FPC 

In this section we present our main technical result, which is an FPC reduction 
from optimisation to separation, which treats the classical polynomial-time re- 
duction of the corresponding problems as a subroutine. This classical result can 
be stated as follows. 

Theorem 8 (c.f., e.g., |GLS88[ Theorem 6.4.9 rj). The linear optimisation prob- 
lem can be solved in polynomial time for any well- described polytope given by a 
polynomial-time oracle solving the separation problem for that polytope. 

Below, we prove the following analogous result for fixed-point logic with count- 
ing. 

Theorem 9 (Optimisation to Separation). Let V be a class of well-described 
rational polytopes represented by T-structures and the function v . Let S be an 
FPC interpretation of r^ec in t ^ r^ec expressing the separation problem for V 
with respect to v. Then there is an FPC interpretation of T(j tt) T„ec m r l+l r, 
which expresses the linear optimisation problem for V with respect to v. 



vec 



Observe that these theorems do not imply that every linear optimisation 
problem can be solved in FPC (or even in polynomial time). Rather one can 
solve particular classes of linear optimisation problems where domain knowledge 
can be used to solve the separation problem. We have the following generic con- 
sequence in the case of explicitly-given polytopes when Theorem |9] is combined 
with Theorem [T] 

Theorem 10 (Explicit Optimisation). There is an FPC -interpretation o/tqI+) 
Tyec in Tmat^Tvec^Tvec cxprcssing the linear optimisation problem for the class of 
polytopes explicitly given by constraints and represented naturally as Tmat^Tyec- 
structures. 

The main idea behind the proof of Theorem |9] is as follows. Suppose we 
are given a polytope P C Q^ by an FPC-interpretation Sp that expresses the 
separation problem for P. A priori the elements of V are indistinguishable. 
However, Ep may expose an underlying order in V as it expresses answers to 
the separation problem for P. For example, suppose Sp on some input expresses 
a vector d € Q^ where the components d„ and d^ for u,v (1 V are different. 
This information can be used to distinguish the components u and v; moreover, 
it can be used to order the components because d^ and d^ are distinct elements 
of a field with a total order. As Sp is repeatedly used it may expose more and 
more information about the asymmetry of P. This partial information can be 
represented by maintaining a sequence of equivalence classes (14)f=i partitioning 
V . This equivalence relation is progressively refined through further invocations 
of the separation oracle. Initially all elements of V reside in a single class. 



^The reverse of this theorem also holds: An oracle for the linear optimisation problem can 
be used to solve the separation problem. 
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It is natural to consider the polytope P' derived from P by taking its quo- 
tient under the equivalence relation defined in this way. Intuitively, this maps 
polytopes in Q^ to polytopes in Q*"' by summing the components in each equiv- 
alence class to form a single new component which is ordered by the sequence. 
We call this process folding. We observe that a separation oracle for P' can 
be constructed using Tip, provided the answers of Ep never expose more asym- 
metry than was used to derive P' . However, failing to meet this proviso is 
informative — it further distinguishes the elements of V — and refines the se- 
quence of equivalences classes. 

These observations suggest the following algorithm. Start with a sequence 
{Vi) of exactly one class which contains all of V. Construct the folded poly- 
tope P' with respect to this sequence and the associated separation oracle Sp' 
from Ep. Attempt to solve the linear optimisation problem on the folded 
polytope (which lies in an ordered space) using the Immerman-Vardi theorem 
[Var82[ IImm86j and the classical polynomial-time reduction from optimisation 
to separation (Theorem plr] Should Ep at any point answer with a vector 
that distinguishes more elements of V than the current sequence of equivalence 
classes, then we: (i) abort the run; (ii) refine the equivalence classes and the 
folded polytope with this new information; and, finally, (iii) restart the opti- 
misation procedure on this more representative problem instance. Since the 
number of equivalence classes increases each time the algorithm aborts, it even- 
tually solves the optimisation problem for some P' without aborting. We argue 
that this solution for P' can be translated into a solution for P. 

A key aspect of this approach is that it treats the polynomial-time reduction 
from optimisation to separation as a blackbox, i.e., it assumes nothing about 
how the reduction works internally|j Before formally describing the algorithm 
we establish a number of useful definitions and technical properties. 

4.1 Folding 

Let y be a set. For k < \V\, let cr : y —> [k] be an onto map. We call a an 
index map. For i £ [k] define Vi := {s £ V \ a{s) = i}. The sequence of sets Vi 
is a partition of V . 

Definition 11 (Folding). 

For a vector x € Q^ let the almost-folded vector [xY ofQ'' be given by 



i[x]'^)i := ^ x^,, fori € [k]. 



For a vector x € Q^ let the folded vector [x]"^ of Q'' be given by 

([x]-),:=|^, /or^e[fc]. 
For a vector X £ Q*^ let the unfolded vector [x]^'^ of Q^ be given by 

(M^°^)f ■— ^ij ™*^^ ^i 3 ^' /""^ V £V. 



^Recall that, by the Immemian-Vardi theorem, every polynomial-time property of ordered 
structures is definable in fixed-point logic, and hence also in FPC. 

*It is, in fact, possible to translate the classical reduction from optimisation to separation 
line by line into FPC (in the spirit of Sectional. However, this translation quickly becomes 
mired in intricate error analysis which is both tedious and opaque. 
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We say a vector x S Q^ agrees with a when for all v,v' G V, (t{v) = (t{v') 
implies Xy — Xyi. It easily follows that if x agrees with a then [[2:]°^]^'^ = x. 
When vectors agree with a and a is clear from context we often use the font, 
as above with x and x, to indicate whether a vector is unfolded and lies in 
Q^, or folded and lies in Q*^, respectively. The notion of folding naturally 
extends to a set 5 C Q^ (and hence polytopes): Let [S*]"" := {[s]"^ | s G 5}. See 
Figures [2] and [3] for examples of folding polytopes. Note that [P]'^ is a projection 
of P into the /c-dimensional space Q*"' . 

Several useful properties of folding and unfolding follow directly from their 
definitions. 

Proposition 12. Let a : V ^ [k] be an index map and c,x E Q such that c 
agrees with a . Then, 

c^Wxrv" ^c^x^\cf'^\xr. 



Proof. We begin by proving the first equality. Fix i e [k]. Definition 11 implies 
that 

J2 mi-^v -^v^J2 ( W')^ - E ^^ = E ( r^ E ^-' ) - E ^- 



WA 



/ ^ ^v' / ^ '^v — U- 



v'eVi veVi 

We conclude that 

c'^{[[xr]-''-x) = j2cM[[<]-i.~^.) 

vev 

= E E ^"(([N^-^). - 

ie[k] veVi 



(1) 



V = ^^e[k]V^ 



E E ^^iiii<rnv - x,) c agrees with a 



iG[fe] vl^Vi 



E ¥?^ E [(MT^). 



ie[k\ 



veVi 



E&-o = o. 



i^[k\ 



VA 



linearity 




We now argue the second equality. 
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Figure 2: Folding and unfolding a polytope P C Q^^"'"!^ with respect to a 
{u-)-0,v-)-0}. 
w 




~ U — V 



Figure 3: Folding and unfolding a polytope P C qI"-^-"'} with respect to a 
{u —i' 0,v —i' 0,w ~> 1}. 
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c^ 



vev 

ie[A;] -ueVi 

= ^ ^ I jTTj ^ c„' I Xy c agrees with ct 

-E(W^)<:(N^). Def.[il] 

■ie[fc] 



n 



4.2 Folding Polytopes 

The diagrams in Figures [2] and |3] suggest intuitively that the result of folding 
a polytope is itself a polytope; the following proposition makes this connection 
concrete. 

Proposition 13. Let P be a polytope in Q^ and let a : V -^ [k] he an index 
map. Then the folded set [P^ is a polytope with {[PY)i < A8k^\V\^{P)f. 

Proof. Let P = conv(S'i) + cone(S'2) for two finite sets of points 5i, 5*2 C Q^. 
By the linearity of [■]"' we have 

[P]" = [conv(5i) + cone(S'2)]'' 
= [conv(5i)]" + [conv(52)]" 
= conv([5i]") + cone([52]"). 

We conclude that [P]'^ is a polytope. We have 

{[PV)f < 3fc'([P]")v < 3A:2 . 4fc|F|(P)v < 12k^\V\ ■ 4|y|2(F)f 

where the middle inequality comes from bounding the bit complexity of -ryj- for 
extremal vertices v € P. D 

For a polytope P C Q^ and a point x e Q^ (with x ^ P) we say that all 
separating hyperplanes at x disagree with a if there is no c G Q^ which both 
agrees with a and has c^ x > maxjc^y | y e P}. This induces an alternative 
characterisation of the polytope [P]'^ . 

Lemma 14. Let P be a polytope in Q^ and a : V ^ [k] be an index map. Then 

\p]c _ p' J (n,k W '^ £ P or all separating hyperplanes 

\ at [x]""^ disagree with a 
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Proof. We show both inclusions. 

1. [P]*^ C P': Let X e [P]". By definition there is a point x & P such that 

[x]°' = X. Suppose X — [x]"'^, then [x]"'^ G P and hence x e P'. Thus assume 



[x]^'^ 7^ a;. Let c G Q^ be any vector agreeing with a. By Proposition 12 we have 
c^[x]^'^ = c^[[a;]'^]^'^ = c^x. Since x G P, c is not the normal of a separating 
hypcrplane through [x]"*^. We conclude that all separating hyperplanes through 
[x]^'^ disagree with a and hence that x E P' . 

2. [P]" 2 P': Let x e P'. Suppose [x]-'" e P, then x = [[x] "'"]'" G [P]'". 
Thus assume that [x]^'^ ^ P and that all separating hyperplanes through [x]^*^ 
disagree with a. This means that for any vector c e Q^ that agrees with a 
the hyperplane through [x]""' with normal c intersects P and thus there is a 
point y E P which has c^[x]^°" = c^y. This further implies that c^[x]^°' < 



max{c y I y G P}. Since c agrees with a, Proposition 12 implies that 



[cf^x<max{[cf^[yr\yeP}. 



Observe [{c G Q^ | c agrees with cr}]"' = Q'^. This means for any vector c' e Q*^, 
c'^x < max{c'^[y]'^ | y £ P}. In particular, for every constraint defining the 
polytope [P]'^, X also satisfies that constraint. We conclude that x G [P]"^. D 

4.3 Expressing Optimisation in FPC 

Suppose we are given a polytope P C Q via a separation oracle Ap, and a 
vector c indicating a linear objective. The algorithm maintains an index map 
a : V —>■ [k] that indicates a sequence of equivalence classes of V which have not 
been distinguished by the algorithm so far. Initially this index map is given by 
ordering variables according to their relative values in c. Under the assumption 
that a accurately describes the symmetries of P we execute the polynomial-time 
reduction from optimisation to separation on the polytope [P]"' and objective 
[c]'^. Since [P]"' lies in an ordered space, it follows from the Immerman-Vardi 
theorem that the reduction can be expressed in fixed-point logic with counting. 
To this end, a separation oracle A[p]<t must be specified for the polytope 
[P]°". Given a point x G Q*^, we argue that the result of applying Ap to the 
unfolding of x either determines the point is in P, and hence also in [P]"'; 
or determines a separating hyperplane for P. If a separating hyperplane is 
determined, it can be folded into a separating hyperplane for [P]'^, but only if 
the hyperplane normal agrees with a. In the case the separating hyperplane 
disagrees with a, our assumption about P is violated, and our separation oracle 
does not have enough information to proceed. Indeed, folding the resulting 
normal may produce 0*^ which is not a valid answer. In this case, the algorithm 
aborts the run of the linear optimisation algorithm, and returns the disagreeing 
hyperplane normal. The algorithm then combines the disagreeing normal with 
its current index map a to produce a new index map which is consistent with a 
and agreed with by the disagreeable hyperplane normal. This strictly increases 
the number of equivalence classes of variables induced by the index map. The 
above procedure can abort at most \V\ times before a exactly characterises the 
order of V relative to P. After this point the linear optimisation algorithm 
cannot abort and hence must solve the optimisation problem for [P]"' which can 
be unfolded into a solution for P. 
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Figure 4: An instninientation of the reduction from optimisation to separation. 
OPT*(P,Ap,c) 
Input: 

• A well-described polytope P C Q^ with a separation oracle Ap, and 

• a linear objective c G Q^. 
Output: 

• / = 1 and y = 0^ , if P is unbounded along c; or otherwise 

• / = and y e Q^ , s.t. if P 7^ then y G P and c^y = maxjc^x | x e P}. 



1: 
2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 
10: 
11: 



a ^— Refine(0^,c). 
while true do 

(/,x)^Opt([P]-,A[p].,[c]'^). 

if aborted with a' then 
a <— cr'. 

else 

cr' ^ REFINE(cr, Ap([x]"'^)). 

if (T 7^ (t' then 

a -(— cr'. 
else 

return (/, [x]^""). 



12: 
13: 
14: 
15: 
16: 
17: 



oracle Ajp]<t(x) 

d^Ap([x]-). 

a' -It- REFlNE(cr, d). 

if CT 7^ cr' then abort with cr'. 

return [d]"' . 
end oracle 



With this intuition in mind the formal proof is as follows. 

Proof of Theorem^ For completeness the entire algorithm Opt* is described in 
Figure|4j The algorithm uses two subroutines Refine and Opt. The subroutine 
Refine (cr, d) takes as input an index map cr of T^ represented in N^ and a vector 
d € Q^ and computes a new index map a' with the following two properties: 

• for all v^v' ^V with a{v) < <j{v'), <j'{v) < (j'{v'), and 

• for all v,v' eV with a{v) = cr(ti'), a'{v) < cf'{v') iff d^ < dy>. 

It is straightforward to observe that when REFiNE(cr, d) produces an index 
map a' which is different from a, then ct' induces strictly more equivalence 
classes on V then a does. Clearly, no index map can induce more than \V\ 
equivalences classes. The subroutine Opt solves the linear optimisation problem 
on an ordered space Q*^ with a given linear objective and a polytope given by 
a separation oracle. Without loss of generality assume Opt returns a integer- 
vector pair {f,y) which is (1,0'^) when the objective value is unbounded and 



17 



(0, y) when y £ Q^ is an optimal point in the polytope if, and only if, the 
polytope is non-empty. 

We first argue that the algorithm is correct, assuming the correctness of 
OPT and Refine. For any index map cr : V -^ [k], [P]°' is a polytope by 

to 
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Proposition 13 We show that the procedure A[p]<T described in lines 

acts as a separation oracle for [P]" provided the answer given by the separation 

oracle Ap agrees with cr. If Ap([x]"'^) outputs d = 0^, then this indicates 



that [x] "^ € P, and hence x e [P]°' by Proposition 14 Trivially 0^ agrees 
with (J, so [d]°' = [0^]'^ = 0*^ is returned by A[p].t correctly indicating that 
X e [P]''. Otherwise, d ^ 0^ and indicates that [x]"'^ ^ P but d^[x]"'^ > 
maxjd^y | y G P}. If d agrees with a we have, by Propositio n |l2[ [d]'^ x > 
ma.x{[d]'^^ [y]" \yeP}. This is equivalent to [df^x > max{[d]* y | y G [P]"}- 
Hence [d]°' is the normal of a separating hyperplane of [P]°' through x. Since 
d agrees with a, a' = a and [d]"" is correctly returned. If d does not agree with 
cr, then Refine produces a cr' ^ cr and the procedure aborts. We conclude that 
(i) when A[p],T does not abort it behaves as a separation oracle for [P]"^ , and 
(ii) when Arpi.^ aborts the returned index map a' is a strict refinement of a. 
Thus A[p]<T is a separation oracle for [P]"' , provided it does not abort. When 
OPT runs on A[p]<t without aborting the result must be a solution to the linear 
optimisation problem on [P]"^ . 

Let X e [P]" be such that [c]*^x > max{[c]*^y | y G [P]"}, i.e., it is a 
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solution to the linear optimisation problem on [P]"' along [c]'^. By Proposition 
this means that either (i) [x]^'^ e P or (ii) Ap([x]^'^) must disagree with 

a. Applying Ap to [x]^'^ distinguishes these two cases. In case (i), [c]"' x = 
c^[x]~'^ > maxjc^y | y £ P} by Proposition 12 because the initialisation of 



a forces c to agree with cr. This means that [x]""' is a solution to the linear 
optimisation problem for the polytope P and the objective c. In case (ii), 
[x\~°' ^ P but Ap([x]~'^) is guaranteed to improve the index map. In the 
case that the linear optimisation algorithm returns that [P]"' is unbounded in 
the direction of [c]'^ , it implies, via similar analysis, that P is unbounded in the 
direction c. Finally, when the optimisation algorithm reports that [P]'^ is empty 
we conclude that P must be empty as well, because if P contains at least one 
point then [P]"' must also contain at least one point. The algorithm correctly 
translates the solutions for the linear optimisation problem for [P]"' back to 
solutions for P. This means that when Opt* returns its result is correct. 

We now observe that this algorithm runs in polynomial time. The main loop 
cannot execute more than \V\ times, because, as established above, at each step 
either the index map a is improved to induce more equivalence classes — up to |T^| 
classes — or the algorithm returns a correct solution to the linear optimisation 
problem on P. The size of all of the objects referred to by the algorithm can 
be polynomially bounded by a function of the input length. In particular, since 
P is well-described by Ap, there is a polynomial bound on its bit complexity 



and this induces a bound on the size of [P]°' through Proposition 13 and implies 
that [P]"' is well described. This implies that the bit complexity of values in 
the algorithm can be bounded by some fixed polynomial. This means that 
folding and unfolding can be computed in polynomial time. Similarly, a naive 
implementation of the subroutine Refine can be seen to run in polynomial 
time in \V\ and the bit complexity of its input rational vector. Since [P]'^ is a 
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well-described polytope with a polynomial-tinie separation oracle Arpi^ we can 
use the polynomial-time algorithm for OPT from Theorem |8] to solve the linear 
optimisation problem on [P]"^ ■ Combining all these parts implies that Opt* is 
a polynomial-time algorithm. 

We conclude by arguing that the behavior of Opt* can be simulated in FPC. 
Relative to an index map a expressible in FPC, folding and unfolding can be 
expressed in FPC using basic rational arithmetic. It is similarly routine to ex- 
press Refine in FPC by defining the equivalence classes and then counting sizes 
to determine the correct position of each equivalence class relative to an FPC- 
definable a and vector. Moreover, there is a FPC-interpretation S[p]^ expressing 
the separation problem for [P]'^ ■ This implies there is an FPC-interpretation 
for the combination of Opt and the separation oracle given by S[p]<^, because 
the polytope [P]"' lies in an ordered space and the Immerman-Vardi theorem 
[Var82[ IImm86| indicates that any polynomial-time property of ordered struc- 
tures can be defined in fixed-point logic (and hence in FPC). It is easy to see 
that the algorithm's main loop and control structure can be simulated in FPC. 
Combining everything gives an FPC-interpretation simulating Opt* and hence 
expressing the linear optimisation problem for P given a FPC-interpretation 
expressing the separation problem for P. D 

In the next four sections we demonstrate a number of applications of Theorem [9] 
for expressing classical combinatorial optimisation problems in fixed-point logic 
with counting. 

5 Application: Maximum Flow 

Let G ~ (y, c) be a graph with non-negative edge capacities, that is, c : VxV ^ 
Q>o. For a pair of distinct vertices s,t € V an (s,t)-/?owis a function / : 
VxV —i' Q>o satisfying capacity constraints < f{u, v) < c{u, v) on each pair of 
distinct u,v G V and conservation constraints X^^ey (■^(^' ^)~/(^' '*')) = on all 
vertices u e V^\{s, t}. The value val(/) of the flow / is simply the difference in in- 
flow and out-flow at t, i.e., J2vevif(^' ^) ~ /(^' ^))- Observe that any flow / can 
be normalised to /' so that for any pair of distinct u,v &V a,t least one of /'(u, v) 
and f'(v, u) is zero (i.e., if f{u, v) > f(v, u), set /'(w, v) :— f(u, v) — f{v, u) and 
f'{v,u) := 0; obviously this preserves the capacity constraints, the conservation 
constraints and the value of the flow). A maximum {s,t)-flow of G is a, flow 
whose value is maximum over all (s,i)-flows. 

Observation 15. Fix G — {V,c) and s,t £ V. Let /i,/2 be two {s,t)-flows in 
G. Fix any a € Q with < a < 1, let f :— a ■ fi + (l — a) ■ f2. Then /' is an 
{s,t)-flow of G and val(/') — a ■ val(/i) + (1 — a) • val(/2). In particular, if fi 
and /2 are maximum (s,t) -flows then so is any convex combination /'. 

Let G\f := {V, c — f) denote the residual graph of G with respect to the flow 
/. The standard formulation of the maximum (s,t)-flow problem as a linear 
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program is as follows: 

max y^ ifiv, t) — f{t, v)) subject to 



vev 

$](/(«,") -/(u,«))=o, yuev\{s,t} 

vev 



(2) 



5.1 Expressing Maximum Flow in FPC 

Observe that there are |y|(|y| — 1) variables in linear program (l2| corresponding 
to /(u, v) for distinct u,v €V . The program has 2\V\'^ — 4 constraints. Both the 
variables and constraints can be indexed by tuples of elements from V. It can 
easily be established that the maximum (s, i)-flow linear program can be defined 
by an FPC interpretation. That is to say, suppose that a capacitated graph 
{V, c) is given as a r,„at-structure with universe V where the rational matrix c e 
Q>g codes the capacities. Then, there is an FPC interpretation from Tmat to 
Tmat tti Tvcc that takcs a capacitated graph (V, c) and a pair s,t € V and explicitly 
expresses a constraint matrix A and vector b encoding the corresponding flow 
polytope. The flow polytope is bounded because each variable is constrained 
from both above and below. Further the flow polytope is nonempty because 
the capacities in G are nonnegative and hence the zero flow is a member of 
the polytope. Thus, because this polytope is explicit. Theorem |10| immediately 
gives an FPC interpretation expressing the optimisation problem on the flow 
polytope. 

Theorem 16. There is an FPC interpretation $(s,i) of Tmat in T„iat which 
takes a t mat- structure coding a capacitated graph G to a t mat- structure coding a 
maximum {s,t)-flow of G. 

Note that as the interpretation $ defines a particular flow, the flow must, 
in some sense, be canonical because it is produced without making any choices. 
Informally, it is a convex combination of maximum flows resulting from the 
consideration of all orderings consistent with the most refined index map de- 
termined by the FPC interpretation of Theorem [9j This is possible because of 
Observation |15| In our remaining applications — minimum cut and maximum 
matching — the analog of Observation [15] does not hold: Convex combinations 
of cuts or matchings are not necessarily cuts or matchings. In the former it is 
still possible to define the notion of a canonical optimum. In the latter case 
it is easy to observe, as noted in the introduction, that defining a canonical 
maximum matching is not possible. 



6 Application: Minimum Cut 

An {s,t)-cutoi a capacitated graph G — {V,c) is a subset G of the vertices V 
which contains s but not t. The value val(C) of the cut G is the sum of the 
capacity of edges going from vertices in G to vertices in V\G. A minimum (s, t)- 
cut of G is a cut whose value is the minimum over all (s,i)-cuts. A minimum, 
cut of G is a, minimum (s,f)-cut over all choices of distinct vertices s,t. By 
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the max-flow/min-cut theorem, a maximum (s,i)-flow and a minimum. (s,i)- 
cut have the same value. This duahty allows the construction of minimum cuts 
from maximum flows. In this section we describe an FPC formula defining a 
minimum (s, t)-cut in a graph using the FPC interpretation for the maximum 



(s, f )-flow problem given by Theorem 16 we show this minimum cut is canonical 
in a strong sense. 

6.1 Expressing Canonical Minimum Cut in FPC 

First, we define a notion of directed reachability in capacitated graphs. A 
vertex v is reachable from a vertex u if there is a path in the graph which 
follows directed edges with non-zero capacity (this is exactly directed reach- 
ability in the graph induced by eliminating zero capacity edges). Let / be 
a maximum (s,i)-flow in G = (1^, c) with normalised flow /'. Define Cf :— 
{v £ V \ V reachable from s in G|/'}. Cf is a minimum (s, t)-cut in G. Since /' 
is normalised, every edge leaving Cf must be at full capacity in /'. 



Given the FPC interpretation $ from Theorem 16 expressing an (s,t)-flow 
/, it is not difficult to construct a formula of FPC which defines the normalised 
flow /' and then the set of vertices Cf . 

Theorem 17. There is a formula ^(a;, s, t) of FPC which given a Tj^at- structure 
coding a capacitated graph C = (V,c), defines the vertices in a minimum (s,i)- 
cut of G. 

In fact, the cut Cf does not depend on / at all, as wc show next. Indeed, Cf 
is the smallest minimum (s,t)-cut in the sense that it is contained in all other 
minimum (s,i)-cuts of C. 

Lemma 18. Let G = {V, c) he a capacitated graph with distinct vertices s,t £V . 
Then the cut Cf is independent of the choice of a maximum {s,t)-flow f of G. 
Moreover, Cf is the intersection of all minimum {s,t)-cuts of G. 

Proof of Lemma \18\ Suppose not. There are two distinct minimum (s,t)-cuts 
C := Cf and C := Cf with corresponding normalised (s,i)-flows / and /'. 
Since C and C" are different there exists, without loss of generality, v e C'\C. 
Consider the flows through C n C . We use a,a' ,b,b' ,c,c' to denote the net 
flows into and out of this set. See Figure [5] for definitions. By definition of C 
and C there is no flow in / from C to C nor is there flow in /' from C to C" 
as otherwise vertices in the complementary cuts would be reachable from s. 

The flow conservation constraints require the flow into C n C be matched 
by the outflow in both / and /'. This implies that a + b = c and a' = b' + c' . In 
addition a> a' and c' > c, because these edges must be at full capacity in / and 
/' respectively. Combining these equalities and inequalities produces a-\-b < c' 
and b' +c' < a. Adding these two constraints together gives a+b+b' + c' < a+c' . 
Since all values are non-negative we have b = b' — 0. This implies a = c and 
a' — c' . Then, reusing a > a' and c' > c we conclude a — c = a' — c' . This 
means that in /' the edges going from C n C" to C n C" are at full capacity, and 
thus no vertex in G n C is reachable from s in flow /'. This is a contradiction. 

Since the flow from between C D C and G n G' is the same in both / and 
/', flow /' witnesses that G n G' is a minimum (s, i)-cut of G. This implies the 
"moreover" part of the statement and completes the proof. D 
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Figure 5: Diagram for the proof of Lemma [TSJ Here the variables indicate the 
net flow between two sets under flows / and /'. 



Note that this proof is similar to the "lemma on a quadrangle" from |DKL76| , 
but that proof does not immediately go through because C n C" may not be a 
(s,i)-cut. 

Lemma [T8| implies that the FPC formula ^ of Theorem [T7| defines a unique 
(s,t)-cut of the graph G and for this reason we call it the canonical minimum 
(s,t)-cut of G: KG.s,t- 



7 Application: Minimum Odd Cut 

The minimum odd cut problem is closely related to the minimum (s, i)-cut 
problem. Here the goal is to define a minimum odd cut of a graph G. That 
is, a cut of odd size whose value is minimum among all odd size cuts of G. A 
capacitated graph G = {F,c) is sym,m,etric if for all u,v G V, c{u,v) — c{v,u). 
Observe that in symmetric graphs the set C is a minimum (s, i)-cut iff its 
complement C := V\C is a minimum (i, s)-cut. In this section we prove that in 
each symmetric graph G there is at least one pair of vertices s, t such that the 
canonical minimum (s,i)-cut KG^s,t is a minimum odd cut of G. The results 
and techniques discussed in this section are entirely graph theoretic. 

Before continuing we must define several special types of cuts. We extend 
a capacitated graph G = {V, c) to a marked capacitated graph G" = {V, c, M) 
with a marking M C y. We call a vertex v G V marked if w € M . A cut 
G of a marked graph G is said to be a marked cut, if both G and G contain 
a marked vertex. A cut G of a graph G = {V,c) with \V\ even is said to be 
an odd cut if \C\ is odd. A marked cut G of a marked graph G — [V,c,M) 
with \M\ even is said to be an odd marked cut if |G n M\ is odd (note that 
this corresponds to the simpler notion when M = V). For any set C of cuts we 
define the basic cuts in C to be {G e C | VG' e C, G = G' or G ^ G'}. Note that 
if C is non-empty it must contain at least one basic cut. When C is the set of 



minimum (s,i)-cuts, the formula ^ of Theorem 17 defines the unique basic cut 



in C. We frequently describe sets of cuts by a sequence of the above adjectives 
and determine meaning by first evaluating the adjective which appear closest to 
the word "cut" . The most complex cuts we consider are "basic minimum odd 
marked cuts" . 



22 







c 


a 








C" 


h 




i 

• 






a e 
\ 


f c 

/ 


C 






/ 

d 


\ 




c 










•s 

















Figure 6: Diagram for the proof of Lemma [19] 



Section 7.1 develops several technical properties of cuts of marked symmet- 



ric graphs. Using these properties Section 7.2 shows that there is a canonical 
minimum (s,i)-cut of a symmetric graph G which is also a minimum odd cut 
of G. 



7.1 Intersections of Minimum Cuts 

We now prove two technical properties involving the intersections of marked 
cuts. The first says that basic minimum marked cuts do not have complicated 
intersections with basic minimum (s,f)-cuts. The proof is similar in spirit to 
the proof of Lemma [T8| 

Lemma 19. Let G — {V,c^M) he a marked symmetric graph. Let s,t £ AI 
he distinct vertices and let C he a hasic minimum {s,t)-cut of G. For every 
hasic minimum marked cut C" of G one of the following holds: (i) C 3 C , (ii) 
CnC = % or (Hi) {s, t} n C" 7^ 0. 

Proof. Fix any basic minimum marked cut C" of G. Suppose neither property 
(ii) or (iii) holds; it suffices to show that property (i) holds. Thus our goal is 
to show that C" := C n C" = assuming that C n C" 7^ and {s, t)C^C' = 0. 
See Figure |6] for a diagram of the general configuration of these cuts and for the 
definitions of variables labelling the symmetric capacity crossing between the 
various sets. 

Observe that Cn C is an (s, i)-cut, however, it cannot be a minimum (s, t)- 
cut because C 3 Cn C (since Cn C" ^0) and C is a basic minimum (s, i)-cut. 
Thus _ 

c + d + e = val(C n C) > val(C) = a + c + e + f 

and hence d > a + f. 

Since C is a marked cut of G, C contains a marked vertex. Suppose that 
C" contains a marked vertex. In this case C" is marked cut of G, because it 
contains at least one marked vertex, but not all marked vertices (e.g., s). Since 
C is a basic minimum marked cut of G and C" C C', C" cannot be a minimum 
marked cut of G. This implies that 

a + 6 + e = val(C") > val(C') =b + d + e + f 
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and hence that a > d + f. Combining this with the inequality d > a + f derived 
from C being a basic minimum (s,t)-cut we have d>a + f>d-\-2f which 
is a contradiction because all the edges have non-negative capacity. Thus C" 
cannot contain a marked vertex. This implies that C n C" contains a marked 
vertex and is hence a marked cut of G. 

Suppose that C" ^ 0. Since C" is a basic minimum marked cut of G, 
C n C" C C" cannot be a minimum marked cut of G. Thus 

a + d + f = val(C n C") > val(C") = h + d + e + f 

and hence a > 6 + e. Combining this with the value of G and the fact that the 
edges are non-negative implies 

val(C) =a + c + e + f>b + c + 2e + J>b + c + f^ val(C n C'). 

Thus val(C) > val(C' n C'). As (7 n C' contains t but not s, it is an (i, s)- 
cut. Moreover, C n C is an (i, s)-cut with value strictly less than that of G 
which is a contradiction because C is a minimum (s,t)-cut and the capacities 
are symmetric. Therefore G" — and the proof is complete. D 

The second lemma says that given a basic minimum odd marked cut C, there 
exists a minimum marked cut C" which does not have a complicated intersection 
with G . The proof is quite similar to those of Lemmas [T8| fc [T9l 

Lemma 20. Let G = (V, c, M) be a marked symmetric graph with \M\ even. Let 
G be a basic minimum odd marked cut of G. There exists a minimum marked 
cut G' of G such that one of the following holds: (i) G ^ G' or (ii) G r\G' = ^. 

Proof. Fix any minimum marked cut C" of G. The complementary cut C of 
C" is also a minimum marked cut of G because G is symmetric. If C 3 C" or 
C n C" = the condition is immediately satisfied. If C" 3 C, then C has no 
intersection with G and the minimum marked cut C satisfies (ii). If CnC = 0, 
then G ^ G' and hence minimum marked cut C satisfies (i). In the case that 
none of these things happen we observe that CnC", CnC, CnC and CnC 
are all non-empty. Because C is an odd marked cut there are two disjoint (but 
symmetric) cases: 

1. CnC is an odd marked cut. 

See Figure [7] for a diagram of this case and for the definitions of variables 
labelling the edges crossing between the various sets. The are two further sub- 
cases: 

l.a. CnC contains a marked vertex. 

Since C is a marked cut of G, C contains a marked vertex. Since CnC 
also contains an marked vertex, C n C is a marked cut of G. Furthermore, C 
is a minimum marked cut of C and hence 

b + d + e + f = val(C) < val(C n C) = 6 + c + /. 

We resolve that c> d + e. Similarly, since C is basic, the odd cut C n C C C 
of G must have a larger value than C, and thus d > c + e. Combining the two 
inequalities we conclude d>c + e>d + 2e. This is a contradiction because 
values are non-negative. 
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Figure 7: Diagram for the proof of Lemma |20[ 



l.b. C n C contains no marked vertices. 

Since C is a marked cut, it contains a marked vertex. This imphes that 
C" := C n C contains a marked vertex, and hence that C" is a marked cut of 
G. Thus to satisfy (i) it suffices to show that C" is a minimum marked cut of 
G. 

Because C is an odd marked cut and C n C" is an odd marked cut, C" 
contains an even number of marked vertices. Since C n C contains no marked 
vertices, C contains an even number of marked vertices. This imphes that 
C n C" is an odd marked cut of G. Since C is a minimum odd marked cut of G 
and C n C" is an odd marked cut we have that 

a + c + e + f ^ val(C) < val(C n C") = a + fo + e. 

This implies that h > c + f. Similarly, since C" is a minimum marked cut of 
G and G" is a marked cut of G we have that c > b + f. Combining these two 
inequalities b > c+f > b-\-e+f. As all values are non-negative we must conclude 
that b = c and e = / = 0. This means that the value of G" is 

val(C") = c + d + f = b + d + e + f = val(C") 

and we conclude that C" is a minimum marked cut of G. 
2. G C\G' is an odd m,arked cut. 

Repeat Case 1 with G' and C swapped. D 

7.2 Some Canonical Min (s,t)-Cuts are Min Odd Cuts 

In this subsection we show that for any symmetric graph G ~ (V, c) there exists 
a pair of vertices s,t € y such that the canonical minimum (s,t)-cut Kq^s.i is 
a minimum odd cut of G. 

The intuition for the proof is as follows. Let G be a symmetric graph with 
minimum odd cut G . Mark all vertices in G. We use Lemma [20] to locate a 
basic minimum marked cut D oi G which G does not partition (that is, D is 
either contained in G or disjoint from C). 11 D is an odd cut, and because D is 
basic, a canonical minimum (s, i)-cut separating a vertex s £ D from a vertex 
t ^ D is also a minimum odd cut of G and we are done. Otherwise |-D| is even 
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Figure 8: A subroutine to collapse a set of vertices in a graph to a new single 

vertex. 

Collapse(G, £>, z) 

Input: A marked capacitated graph G ~ {V, c, M), a set D C V and z ^ V. 
Output: The graph obtained from G by collapsing of the vertices in D to z. 

1: V ^V{z/D). 

{c{u,w), u,w£V\D, 

J2xeD^i^^^)' ueV\D,w = z, 
Exe_Dc(2^,w), w £V\D,u = z. 
3: M' ^ M\D. 
4: return (V^',c', Af). 



and we form a new graph G" by collapsing D into a new super- vertex z, and 
then setting the effected capacities so that value of cuts which do not partition 
D are unchanged. Since C does not partition D and \D\ is even, the collapsed 
version C" of C is a minimum odd marked cut of G". We repeat this collapsing 
procedure maintaining a graph G' and minimum odd marked cut C" until we 
locate a minimum marked cut D of G' that is also minimum odd marked cut. 
This provides marked vertices s,t such that the canonical minimum (s,i)-cut 
of G' is a minimum odd cut of G' . The proof concludes using Lemma [l9] to 
translate this fact back to the original graph G and in doing so argues that the 
canonical minimum (s,i)-cut of G is a minimum odd cut. 

We now formalise the notion of collapsing a graph. We begin by establishing 
notation for substituting sets into sets. Let C,D CV and z ^ V such that either 
G 3 L» or G n L» = 0. Define C{z/D) to be a subset of V := iV\D) U {z} 

G(z/i.):J(^\^)^^^>' ^^^' 

^ ' ^ \c, cnD = (l), 

and for G' C V' define C'{D/z) to be a subset of V 

^,(^/,)^ f(C'\{z})uA zea, 
^ ' ' \C', z(^C'. 

Observe that {C{z/D)){D/z) = G. 

The subroutine Collapse(G, D, z) in Fig. [s] describes a method of collapsing 
the vertex set D CV in the marked graph G = {V, c, M) to a single new super- 
vertex z. This subroutine is designed to preserve a basic minimum odd marked 
cut G with respect to a well-chosen marked cut D. 

Lemma 21. Let G = (V,c,M) be a marked symmetric graph with \M\ even. 
Let z <^V and G,D CV be marked cuts of G such that ODD or G Ci D — 9- 
Define C' :^G{z/D) and G' := {V',c',M') := Collapse(G, D, z). 

L The value of G in G is identical to the value of C' in G' . 

2. If G is a basic minimum odd marked cut of G and |_D n A/| is even, then 
G' is a basic minimum odd marked cut of G' , and \M'\ is even, non-zero 
and M' C M. 
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Proof. Suppose C ^ D. The definition of Collapse implies that 

val(C) = ^ c{u,v)^ ^ c{u,v)+ ^ c{u,v) 

u€C,veC u£C\D,v£C u£D,v£C 

uec\D,vec vec 

X c(u,w)=val(C"). 

ue{c\D)u{z},vec 

The case of C n Z) = is analogous. We conclude property 1 holds. 

Assume the hypothesis of property 2. Observe that D contains an even 
number of the marked vertices M but not all of them because Z? is a marked 
cut of G. The subroutine Collapse sets M' = M\D. Therefore \M'\ is even 
because |-DnM| and \M\ are even, and thus M' meets the required conditions. 

We also determine that C is an odd marked cut of G', because |D n Af | is 
even and C is an odd marked cut which either contains D or is disjoint from D. 

Suppose C" is not a minimum odd marked cut of G". Then there exists an 
odd marked cut C" of G" with smaller value than C". The cut C"{D / z) is a 
marked cut of G which does not partition D. By property 1 val(C"(£'/z)) = 
val(C") and val(C) = y&\{C{z/D)) = val(C"). Hence 

val(C"(i:>/z)) = val(C") < val(C") = val(C(z/i:>)) ^ val(C). 

However, C"{D/z) is an odd marked cut because C" is an odd marked cut, z 
is not marked and D contains an even number of marked vertices. This means 
that C is not a minimum odd marked cut of G which contradicts the hypothesis. 
Therefore C" is a minim,um. odd marked cut of G". 

Similarly, suppose C" is not a basic minimum odd marked cut of G". Then 
there exists a minimum odd marked cut G" of G' with G" C G' . By property 1, 
it follows that G"{D/ z) is a minimum odd marked cut of G with G"{D/ z) C G, 
this contradicts the basicness of G . Therefore G' is a basic minimum odd marked 
cut of G' and the proof is complete. 

n 

With the key properties of Collapse established we are ready to prove the 
main result of this section. 

Theorem 22. Let G — (V, c, M) be a m,arked symmetric graph with even \M\ > 
0. There exist s,t £ M such that the canonical minimum {s,t)-cut of G is a 
minimum odd marked cut of G. 

Proof. Fix a basic minimum odd marked cut G of G. Figure |9] describes the 
algorithm WitMinOddCut that computes vertices s and t witnessing the claim 
of the theorem from G by iteratively collapsing G. To prove the theorem it 
suffices to argue the algorithm halts and produces {s,t) £ V^ such that the 
canonical minimum (s,i)-cut KQ^s,t is a minimum odd marked cut. 

As the algorithm runs it maintains the invariant that G* is a basic minimum 
odd marked cut of G*. Observe that this is initially true for G'^ — G because 
M° = M and G is a basic minimum odd marked cut of G. Suppose G' is a basic 
minimum odd marked cut of G' . The basic minimum marked cut C of G' with 
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Figure 9: An algorithm producing a witness for a minimum odd marked cut. 

WitMinOddCut(G, C) 

Input: A symmetric graph G — {V,c,M) and a minimum odd marked cut C 
oiG. 
Output: {s,t) G M^ such that Kq ^t is a minimum odd marked cut. 

qo :=(i/0,c",MO)^ {V,c,M). 

while true do 

Let D^ be a basic minimum marked cut of G" such that C" 15 D^ or 
C n D' = 0. 

if £)' is an odd marked cut of G" then 

return {s,t) with s e D' n Af* and t e KP\D\ 

G*+i ^ COLLAPSE(G',i:'%z')- 

c*+i ^ c*(zVi:>*). 

i -;— i + 1. 



9: 
10: 



C* 3 D^ or C* n C = is guaranteed to exist by Lemma 20 (if the cut given 
by that lemma is not basic there must be a basic minimum marked cut strictly 
within it that continues to satisfy the intersection properties with C). If D^ is 
an odd marked cut, the algorithm halts at linej?] Otherwise the graph G* and 



cut G' are collapsed relative to D\ The second property of Lemma 21 implies 
that G*^^ is a basic minimum odd marked cut of G*^^. Thus the invariant 
holds. 



Lemma 21 also implies that as the algorithm runs, |M*| is even and M*+^ C 
M\ The invariant and Af° — V imply that the test in line |6] will be successful 
and cause the algorithm to halt within y iterations. Because D* is a marked 
cut of G% when line ^ is reached D^ n M' and AP\D'^ are non-empty disjoint 
sets. This means that distinct s and t exist and are returned by the algorithm. 
Let r be the value of i when the algorithm halts. Fix any s e D^ n M^ and 
t G D^\M^. We use the shorthand K^ to denote the canonical minimum {s,t)- 
cut K^i gi for < z < r. It remains to argue that K^ — KG,s,t is a minimum 
odd marked cut. 

Since the cut D^ is a minimum marked cut of G^ and s,t G M^, D^ is also a 
minimum (s, i)-cut of G^ and it has the same value as the canonical minimum 
(s,i)-cut K^. This implies that K^ is a minimum marked cut of G^ because 



s,t & iW. By Lemma 18 D^ contains K^ , but D^ is also basic, so we conclude 
that D'' = K\ 



The first property of Lemma 21 implies that val(G'(zY-D*)) = val(G*+^) 
for all < i < r, and hence that val(G*) = val(G) for all < « < r. 
Since K^ is a minimum marked cut of G"" and G' is a marked cut of G'', 
Ya\{K'^) < val(G'') = val(G). The first part of Lemma [2l| also implies that 
val{K'+^{Dyz')) = va\{K'+^) for all < « < r. Since K'+\D'/z') is an (s,i)- 
cut of G^ val{K') < valiK'+'^iD'/z^)) for all < i < r. Hence we conclude 
that val(ii:0) < val(i^'~) < val(G). 

It remains to argue that K^ is an odd marked cut. Since K^ is an odd 
marked cut, it suffices to show that K^ is an odd marked cut if i^*+^ is an odd 
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marked cut, for all < i < r. To this end assume that K"^^^ is an odd marked 



cut. Apply Lemma [19| with G\ K' and D'; we note that (i) K^ ^ D\ (ii) 
K^ n £>* = 0, or (iii) {s,t} n -D' ^ 0. Property (iii) cannot hold because s and 
t are selected after Z?* was collapsed. This means that K"^ either contains all of 
D^ or is disjoint from D^ . Because the K^ and _?i'*+^ are canonical (s,i)-cuts, 

K' C K'+^{D'/z') C {K\z' / D')){D' / z') = K\ 

As the algorithm did not halt at step i, jZ?* n lvP\ is even and thus K"^ is an odd 
marked (s,i)-cut. 

We conclude that K^ is an odd marked cut with value at most that of a 
minimum odd marked cut C of G. Therefore K'^ = K^s.t is a minimum odd 
marked cut. D 

Consider the following procedure for locating a set of minimum odd marked 
cuts in a marked symmetric graph G — {V,c,M): For all distinct s,t G M 
compute the canonical minimum (s,i)-cut Kq^^i, eliminate those cuts which 



are not odd, then eliminate those cuts which are not minimal. Theorem 22 
indicates that some cuts remain and that those cuts are minimum odd cuts of 
G. Note that the algorithm WitMinOddCut in the proof of Theorem [22] is 
used only in the analysis and not actually run during the above procedure. This 
simple procedure for defining a non-empty set of minimum odd cuts is critical 
to expressing the separation problem for the matching polytope in FPC. 

8 Application: Maximum Matching 

Let G = (V, E) be an undirected graph. A matching M C E is defined by the 
property that no two edges in M are incident to the same vertex. A matching 
M is maximum if no matchings with size larger than M exist. A maximum 
matching is perfect if every vertex in G is incident to some edge in the matching 
(i.e., \M\ = ^). 

8.1 Maximum Matching Program 

Maximum matching has an elegant representation as a linear program. In fact, 
it is an instance of a slightly more general problem: ^-matching. Let c € Q>oj 
b G N^ and A e {0, ij^^^ be the incidence matrix of the undirected graph 
G = {V,E): the columns of A correspond to the edges E and the rows to the 
vertices V, and A^e = 1 if edge e is incident on vertex v. Alternatively we view 
edges e G E as two-element subsets of V. The goal of the 6-matching problem 
is to determine an optimum of the following integer linear program 



„T 



max c 



y subject to Ay < b, y > 0^ . (3) 



We obtain the usual maximum matching problem in the special case where 
6 = 1^ and c = 1^. 

Gcnerically, integer programming is NP-complete, so instead of trying to di- 
rectly solve the above program we consider the following relaxation as a rational 
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linear program. 

max c y subject to 

Ay < b, 
y > 0^ (4) 

yiw) < hb{W) - 1), \fWCV with b{W) odd, 

where y{W) := Y.eeE,e<zwye and b{W) := Y.vew'^v Here we have added a 
new set of constraints over subsets of the vertices. The integral points which 
satisfy (|3|, also satisfy the additional constraints that are added in Q. To see 
this, let y be a feasible integral solution, consider some set W with b{W) odd. 
If \W\ = 1, then y{W) — because no edges have both endpoints in W, so 
assume \W\ > 2. It follows that 2y{W) < b{W), by summing the constraints of 
Ay < b over W with respect to only the edges with both endpoints in W. Since 
b{W) is odd, ^b{W) is half integral, but y{W) is integral because y is an integral 
solution; this means the constraint y{W) < ^{b{W) — 1) is a valid constraint 
for all integral solutions. In fact [lMm65j shows something stronger. 

Lemma 23 ( |Edm65l Theorem P]). The extremal points of the linear program 
(HI) are integral and are the extremal solutions to the b-matching problem. 

Thus to solve 6-matching it suffices to solve the relaxed linear program (|4|. 
As mentioned before, it will not be possible to show that FPC can generally 
define a particular maximum matching, there can be simply too many. However, 
the above lemma means that the existence of a (likely non-integral) feasible 
point y of (HI) with value c^y witnesses the existence of a maximum 6-matching 
with value at least c^y. In addition, the number of constraints in this linear 
program is exponential in the size of the graph G. Thus, we cannot hope to 
interpret this linear program directly in G, using FPC. Rather what we can 
show is that there is an FPC interpretation which, given G, b and c, expresses 
the separation problem for the 5-matching polytope in the linear program (|4|. 
Combining this with Theorem [9] gives an FPC interpretation expressing the 
6-matching optimum. 

8.2 Expressing Maximum Matching in FPC 

The ^-matching polytopes have a natural representation over Tmatch '■= '''mat W 
Tvec- Although the number of constraints in the 6-matching polytope may be 
large, the individual constraints have size at most a polynomial in the size of 
the matching instance. Thus this representation is well-described. 

We now describe an FPC interpretation expressing the separation problem 
for the 6-matching polytope given a Tmatch-structure coding the matrix A and 
bound vector b. As in the explicit constraint setting, our approach is to come 
up with a definable set of violated constraints iff the candidate point is infeasi- 
ble. We then define a canonical violated constraint by summing this definable 
violated set. Identifying violated vertex and edge constraints can easily be done 
in FPC as before. However, it is not immediately clear how to do this for the 
odd set constraints. 

To overcome this hurdle we follow the approach of |PR82| . Let y be point 
which we wish to separate from the matching polytope. Define s := b — Ay to 
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be the slack in the constraints Ay < b. Analogous to b{W), define s{W) :— 
'Evew^v- Observe that 2y{W) + yiW : V\W) + siW) = b{W) (here y{W : 
y\W^) is sum of edge variables with one endpoint in W and one in V^\Vl^). This 
translates the constraints y{W) < ^{b{W)-l) exactly to y{W : V\W) + s{W) > 
1. This means to find a violated constraint of this type it suffices to find W 
such that y{W : V\W) + s{W) < 1. 

Define a marked symmetric graph H over vertex set U :— VU{z} where z is a 
new vertex. Let H have symmetric capacity d: d{u, v) := j/e when u,v € V and 
M, f e e, and d{u, v) :— Sy when u = z and v £ V. Let M := {v € V \ by is odd}. 
If \M\ is odd, add z to M. Thus we have a marked symmetric graph H — 
{U,d,M). Consider any odd marked cut W of H, without loss of generality 
z ^ W (otherwise, take the complement). Observe that the value of edges 
crossing the cut is exactly y{W : t^\Ty) + s{W); also note that s{W) is odd. 
Thus there is a minimum odd marked cut W oi H with value less than 1 iff 
there is a violated odd set constraint in Q. 

By Theorem [22] there is a violated odd set constraint iff for some s,t e M 
the canonical minimum (s,i)-cut is an minimum odd marked cut with value 



less than 1. We conclude, using Theorem 17 and Lemma 18, that we can 
define a family of violated set constraints within FPC. Summing these defined 
violated constraints produces a canonical violated constraint which must be non- 
trivial by Proposition [6] Thus, as in Theorem [7] there is an FPC interpretation 
expressing the separation problem for the polytope in the linear program Bl. 

Lemma 24. There is an FPC interpretation of Ty^c in Tmatch W Tyec expressing 
the separation problem for the b-matching polytopes with respect to their natural 
representation as Thatch- structures. 

Like the maximum flow problem in Section [Sj the 6-matching polytope is 
both compact and nonempty. By combining Lemma 24 and Theorem l9] with 



respect to the natural well-described representation of 6-matching polytopes, 
we conclude that there is an FPC interpretation expressing the value of the 
maximum 6-matching of a graph. 

Theorem 25. There is an FPC interpretation of tq in Tmatch W Ty^c which 
takes a Tmatch tt* Tyec- structure coding a b-matching polytope P and a vector c to 
a rational number m indicating the value of the maximum b-matching of P with 
respect to c. 

9 Conclusion 

Our main result is that the linear programming problem can be expressed in 
fixed-point logic with counting — indeed, that the linear optimisation problem 
can be expressed in FPC for any class of polytopes for which the separation 
problem can be defined in FPC. As a consequence, we solve an open problem 
of |BGS99j concluding that there is a formula of fixed-point logic with counting 
which defines the size of a maximum 6-matching in a graph. This is one demon- 
stration of the power of the ellipsoid method and linear optimisation that can 
be brought to bear even in the setting of logical definability. From here, there 
are number of natural research directions to consider. 
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Convex programming. A polytope is an instance of much more general 
geometric object: a convex set. The robust nature of the elUpsoid method 
means it has been extended to help solve more general optimisation problems, 
e.g., semi-definite programs and quadratic programs. It seems likely that our 
methods can be extended to these settings. 

Completeness. Linear programming is complete for polynomial time under 
logspace reductions |DLR79) . It follows from our results that it cannot be 
complete for P under logical reductions such as first-order interpretations, since 
this would imply that P is contained in FPC. Could it still be the case that linear 
programming is complete for FPC under such weak reductions? Or perhaps 
FOC reductions? Even if linear programming is not complete, there may be 
other interesting combinatorial problems that can be expressed in FPC via 
reduction to linear programming. There has also been some work examining 
generalisations and improvements to the 5-matching approach we followed (e.g., 
[CF96J ). and it is possible these results can also be replicated in FPC. 

LP hierarchies and integrality gaps. Another intriguing connection be- 
tween counting logics and linear programming is established in |AM121 IG012] 
where it is shown that the hierarchy of Sherali- Adams relaxations |SA90j of the 
graph isomorphism integer program interleaves with equivalence in fc-variable 
logic with counting (C*^). It is suggested |AM12j that inexpressibility results 
for C'^ could be used to derive integrality gaps for such relaxations. It is a 
consequence of the results in this paper that the Sherali- Adams approximations 
of not only isomorphism, but of other combinatorial problems can be expressed 
in FPC. Do our results provide another route to using inexpressibility in FPC 
to prove integrality gaps? 
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